Lift and down-gradient shear-induced diffusion in Red Blood Cell suspensions 
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The distribution of Red Blood Cells in a confined channel flow is inhomogeneous and shows 
a marked depletion near the walls due to a competition between migration away from the walls 
and shear-induced diffusion resulting from interactions between particles. We investigated the lift of 
RBCs in a shear flow near a wall and measured a significant lift velocity despite the tumbling motion 
of cells. We also provide values for the collective and anisotropic shear-induced diffusion of a cloud 
of RBCs, both in the direction of shear and in the direction of vorticity. A generic down-gradient 
subdiffusion characterized by an exponent 1/3 is highlighted. 



PACS numbers: 47.63.-b,47.57.E-,83.50.Xa,83.80.Lz 

Blood is a dense suspension of deformable cells, mainly 
red blood cells (RBCs), making it a complex fluid from a 
rheological viewpoint, and leading to complex flow pat- 
terns in the microcirculation where the diameter of blood 
vessels becomes comparable to cell size. 

In his pioneering work, Poiseuille revealed that blood 
flow in arterioles and venules features a RBC-free plasma 
layer near the vessel wall [TJ [2] . The lubrication effect of 
this depleted layer leads to the Fahraeus-Lindquist effect, 
a decrease of the apparent viscosity of blood in small 
vessels when their diameters become comparable to cell 
size (d < 500 pm) [3] . 

A classic result in low-Reynolds number hydrodynam- 
ics — relevant to blood flow in arterioles and venules [4j|5] 
— is that migration of spherical particles transversally 
to flow direction is prohibited by the linearity and flow- 
reversal symmetry of the Stokes equation [6 j. However, 
the deformability or the non-sphericity of RBCs allow a 
symmetry breaking that may lead to transverse migra- 
tion, be it due to interactions with walls or neighboring 
cells. 

In a shear flow near a wall, lipid vesicles experience 
a lift force that pushes them away [7HTD], at least when 
they are in a tank- treading regime, with steady inclina- 
tion angle. A straightforward question arises: how do 
RBCs, that are usually in a tumbling regime [4] [11] and 
explore all angles, still experience a non-zero average lift 
force ? While many numerical studies have tried to re- 
produce this behavior [T2HT7] , experimental data on this 
basic mechanism are rare [4 or focused on RBCs artifi- 
cially placed in the tank-treading regime [18]. 

This migration of blood cells forms the physical basis 
of the formation of a depleted layer near vessel walls in 
the microcirculation. However this phenomenon alone 
cannot explain the complexity of flow patterns observed 
in the microvasculature, where redistribution processes 
are indeed very frequent since bifurcations are met every 
20 vessel radii [19]. In physiological conditions, blood 



is a very concentrated suspension with a hematocrit up 
to 50%, in which the hydrodynamic interactions between 
cells play a decisive role. The interactions between two 
bodies in flow is a fundamental question in the frame- 
work of suspension dynamics and rheology, even in rather 
dilute suspensions [20j [21]. Unlike smooth and spheri- 
cal particles, rough spheres [22] and deformable particles 
such as drops, bubbles, capsules or vesicles [23-27 , are 
irreversibly shifted after interaction. 

The cumulative effect of these hydrodynamic interac- 
tions is a non-linear and anisotropic shear-induced dif- 
fusion (SID) [22] [24]. The consequences of this SID are 
twofold: repeated collisions of one blood cell with the 
others lead to a random walk which may be important 
for mixing properties of blood flows (self diffusion) , and a 
redistribution of concentration inhomogeneities, that bal- 
ances lift forces (collective or down-gradient diffusion). 
The coefficients characterizing both phenomena are a pri- 
ori different [22]. Investigations of the random walk of 
RBCs in concentrated suspensions [4] [28] [29] provided 
values of the dispersion coefficient two orders of mag- 
nitude higher than Brownian diffusivity. These stud- 
ies were complemented by cell tracking experiments in 
quasi 2D flow [30 . All these studies and most numerical 
works [31] [32] focus on the self diffusivity, however there 
has been no experimental quantification of down-gradient 
diffusion in non- homogeneous suspensions of RBCs and 
whatever the considered particles, experimental charac- 
terizations of down-gradient diffusion are scarce [33] [34] . 

We report on an experimental study on the lift of di- 
luted RBCs in shear flow near a wall and show that tum- 
bling RBCs follow the same scaling laws as tank-treading 
vesicles [3 [10 , capsules [35] and drops [33 . In a different 
experiment, the collective SID of RBCs was investigated, 
providing values of the diffusivity in the vorticity direc- 
tion and in the direction of shear. 

Lift of RBCs in shear flow near a wall - In order 
to avoid screening by sedimentation, measurements were 
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performed in microgravity in CNES and ESA parabolic 
flight campaigns. The procedure and experimental setup 
are detailed in Ref. [10 . We use a Couette shear flow 
chamber with two parallel glass discs, with a gap of 
170 um. The 3D positions of the RBCs which initially lie 
on the bottom disc are captured by digital holographic 
microscopy [36) [37]. 

Blood was collected from healthy donors and washed 
twice in Phosphate Buffer Saline (PBS) and Bovine 
Serum Albumin (BSA). After gentle centrifugation, 
RBCs were dispersed into different fluids: (PBS)+(BSA) 
alone or combined with a mixture of 1% dextran of molec- 
ular weight 1.5 x 10 4 + n% dextran of molecular weight 
2 x 10 6 , n = 3, 4, 5. Corresponding viscosities are respec- 
tively 1.4, 6.1, 9.3 and 13mPa.s (T w 21°C). 

For a neutrally buoyant ellipsoidal lipid vesicle, the 
theoretical drift velocity is given by [7]' z = UjR 3 /z 2 , 
where R is a particle characteristic size, 7 the shear rate, 
and z the distance to the wall. U is a dimensionless drift 
velocity which depends on vesicle shape, and on the inner 
and outer fluid viscosities. It yields the following scaling: 



3UR 3 jt 



(1) 



The evolution of the mean transverse position (z) 3 is 
presented as a function of jt in FigJT] For a given external 
solution all results fall on the same straight line in agree- 
ment with Eq(T] By symmetry, a tumbling rigid object 
should not migrate on average [7J. The non zero lift sug- 
gests that RBC deformability allows symmetry breaking: 
it is stretched when oriented in the direction of the elon- 
gational component of the flow, while it is compressed 
when orthogonal, resulting in an averaged asymmetric 
shape, leading to a migration law similar to the one 
known for a fixed shape and orientation. By increasing 
the external viscosity, the stresses on the RBC membrane 
are higher and lead to increased deformation [11] which 
in turn enhances the lift. RBCs in the 13mPa.s solution 
have viscosity ratio close to 0.3 [38. and are probably very 
close to tank-treading regime [39] , in which case Olla pre- 
dicts a migration with comparable UR 3 = 6.4 um 3 for a 
vesicle of similar (but fixed) shape [7], with long axis 
equal to 7.2 um, the mean diameter of a RBC [40 . As- 
suming x = 72/, we find that in physiological conditions, 
a RBC will migrate by 8 um while travelling 1 cm, a re- 
sult in good agreement with the pioneering result of 4 um 
drift in Poiseuille flow by Goldsmith [4]. 

Shear-induced diffusion in channel flow — The SID 
of a RBC suspension was studied in standard poly- 
dimethylsiloxane microfluidic chips. Thanks to a flow- 
focusing device, a thin layer of RBC suspension is pro- 
duced in a rectangular channel where a buffer (PBS solu- 
tion) flows in the x direction, the gravity direction. The 
RBC cloud is pinched at the entrance in the y direction 
and the direction of observation is z, allowing to record 
the evolution of the RBC cloud in the (x, y) plane (Fig. 
[2|. We restrict the study to moderate shear rate values 
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FIG. 1: (Color online) RBC - wall distance {z) s vs. jt for 
different outer viscosities: (a) 1.4mPa.s; (b) 6.1mPa.s; (c) 
9.3mPa.s; (d) 13mPa.s. Full lines indicate fit to Eq. [I] with 
UR 3 = 0.36, 3.1, 3.2, 5.4 um 3 , respectively. 



bounded by shear rate at the edge 7 max < 340s -1 and 
therefore comparable to physiological ones [4]. 

We first focus on two channels with high aspect ratios: 
width 2dxheight 2ft = 491x53um 2 and 497xl01um 2 . 
Thus, the velocity profile is parabolic across Oz and al- 
most flat in the Oy direction. In this case, due to the 
strong shear in the z direction, the concentration tends 
to homogenize quickly due to diffusion in the plane of 
shear, while diffusion in the vorticity y direction leads to 
the observed widening (Fig. [2|. 

From microscopic images taken with a long exposure 
time, a calibration process based on the Beer-Lambert 
law relates the grey intensity to the local concentration 
profile (j)(x, y). The evolution of the concentration profile 
along x is directly related to the diffusion process through 
the following advection-diffusion equation [41]: 



dx dy 



d_ 

dy \^ dy 



(2) 



where (.) denotes the average over z. 

Here we assume that the concentration of RBCs is ho- 
mogeneous in the z direction. We also suppose that the 
velocity u is the one of an unperturbed Newtonian fluid 
[42] . The diffusivity D = f 3 R 2A f(j) 1S proportional to the 
frequency of pair interactions 70, a straightforward scal- 
ing for shear-induced diffusion due to pair interactions 
[22] [24] . As in Ref. [22] , we denote the dimensionless 
diffusivity in the vorticity direction. 

Rusconi and Stone made similar experiments with 
platelet like particles, with different initial and bound- 
ary conditions [34]. They considered the spreading of 
a concentration step with fixed concentrations at each 
end and found a x 1 ! 2 scaling for the diffusive front. In 
our case, a peak of fixed area spreads and self-similar 



3 




x = x = 30 mm x = 63 mm x = 100 mm £ = 202 mm x = 303mm x = 411mm £ = 598 mm 



FIG. 2: Example of RBC diffusion in a fiat channel. Initial mean volume fraction is around 15%. 



solutions exist under the condition of a widening with a 
x 1 / 3 scaling [41]. The self-similar concentration profile 
is parabolic and one finds the following relation for the 
expected half-width at half-height of the RBC cloud: 



w(x) = w (l + Ax/wl) [ 



(3) 



where the initial peak has width 2wo, A = 2T ^^^ Ar ° and 
Nq = j (j)(x,y)dy is the conserved number of particles 



[IT] . The scaling 



Ax as well as the parabolic 



concentration profiles are nicely recovered in experiments 
for different /i, wq and Nq (Fig. [3|, and the slope A gives 
a direct measurement of fs- Fig. [4ja) shows that for all 
available data in the mid concentration range (<p < 16%) 
A is a linear function of N /h : giving a dimensionless 
diffusivity for RBCs fa = 0.12 ± 0.01, with 2R = 7.2 urn, 
the mean diameter of a RBC. With similar choice for i?, 
fa — 6.9 was found for very flat platelike particles [34] . 
This discrepancy cannot be related to the deformability 
of RBCs: self diffusivity of hardened cells has been shown 
to be of the same order as the one of normal cells [30 . 
However, both discoidal particles are tumbling, thus the 
effective occupied volume is much larger. Replacing <j> by 
<j)V e /V , where V is the particle volume and V e = AttR 3 /3 
this effective volume, we find fa = 0.05 for RBCs and 
fs = 0.18 for platelike particles. These values are now 
comparable. The remaining difference can be attributed 
to the details of the hydrodynamic interactions. For in- 
stance, in Ref. [22] . fs varies from to 0.03 for rough 
spheres with minimal separation going from to 0.0SR. 

At higher volume fractions of RBCs (0 < 30% in 
the experiments), the diffusion should not be a conse- 
quence of pair-wise interactions only, since a RBC inter- 
acts with multiple neighbors, and interactions between 
at least three bodies should also be considered. These 
interactions lead to an additional term proportional to 
(j) 2 in the diffusion coefficient and to a different scaling, 
namely x 1 / 4 if these 3-body interactions were the domi- 
nant effect [41]. However, the noise in the experimental 
data allows a rather good rescaling with exponents be- 
tween 1/3 and 1/4. By forcing a 1/3 exponent, one gets 
an effective diffusivity which increases more than linearly 
with concentration, showing the increasing importance of 
3-body interactions in the diffusive process (Fig. [4^b)). 

An interesting potentiality of this experimental de- 
vice is the possibility to measure both diffusivities fa 
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FIG. 3: (color online) (a): concentration profiles <j>(x,y) in 
four sections of flat channel with 2h = 53 um, for a cloud 
with wq — 41.0 um and iVo — 8.1 um. Full lines show fits 
with parabolic profile, (b) Cloud half-width it; as a function 
of position x along the channel for several initial conditions 
and for two different thicknesses (empty symbols, 2h = 53 um, 
7max = 113 s _1 ; full symbols, 2h = 101 um, 7 max = 211 s _1 ). 
Full lines show linear fit for w 3 as suggested by Eq. [3] 



and fs corresponding to repulsion of interacting cells in 
the plane of shear and in the vorticity direction [24] . 
In a channel with cross section 190x99um 2 , a nearly 
parabolic flow with gradients of velocity in both direc- 
tions y and z was produced. The averaged concentra- 
tion profile observed in the z direction therefore widens 
due to hydrodynamic repulsion in the local shear and 
vorticity directions. By varying the initial position yo 
of the RBC stream, one can vary the weight of the fa 
and fs contributions, with a contribution of fs only for 
yo = and an increasing contribution of fa as the stream 
is moved towards the channel edges. With the addi- 
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tional simplification that in the y direction all particles 
experience the velocities and shear of position y (nar- 
row cloud approximation), one gets an equation similar 
to Eq. [3j with coefficient A that now depends on y^: 

A = ( / 2 2*y + 2-fi/2 ) a ^ i where ui is the partial deriva- 

tive of u according to variable i = y, z at position (yo, 2) 
[4T] . Considering narrow initial clouds (wq — 6um) con- 
sistent with the above simplification, the scaling with ex- 
ponent 1/3 is confirmed by the experiments and the re- 
sulting effective diffusion coefficient A increases with yo 
(Fig. |5j. Within the experimental uncertainties, A does 
not depend on the mean shear rate, though the RBC dy- 
namics and the consequent interaction trajectories might 
be affected by the shear rate value pT[[43j[44] • A fit of the 
data by the expected expression yields fs = 0.07 ± 0.01 
and fi — 1.7 ± 0.1. The fs value is lower than the one 
previously found. Around yo = 0, shear intensity in the 
y direction vanishes, so does collision rate and diffusion 
is expected to be similar to the one observed in the flat 
channel, and controlled by fs. However, wq is finite, 
and the 3D shear also controls the mean orientation of 
RBCs, therefore the detail and intensity of their inter- 
actions, and finally the resulting diffusion coefficient fs 
may be affected. As for drops [24 or rough spheres [22 , 
/2 is found to be larger than / 3 . In the case of drops, 
experiments of Ref. [33 show that fi — 0.2, which is 
comparable to our f^ — 0.77 obtained after rescaling by 
the effective volume. Finally, these down-gradient diffu- 
sion coefficients should be compared to the self diffusion 
coefficients. Self diffusion in the vorticity direction was 
studied in Ref. [30], but the diffusion coefficient in this 
very flat geometry (12 um thick channel) is surprisingly 
found to be independent from concentration. By lack 
of similar scaling, comparison is therefore not possible. 
Self diffusion in the shear direction is characterized by 
D s /j of order 1 um 2 for ~ 40% H El EE]. We 
find f2R 2< P — 9 pm 2 , a consistent result since this down- 
gradient coefficient is expected to be a few times larger: 
6 times for rough spheres [22] and 5 times for drops [33 . 

Conclusion. - Our quantitative investigation of the 
migration of tumbling RBCs in shear flow shows that 
wall-induced lift follows the same scaling law as particles 
with fixed orientation (e.g. tank-treading vesicles), and a 
significant amplitude has been measured. The necessary 
symmetry breaking is made possible by RBC deforma- 
bility. In blood vessels, lift is balanced by shear-induced 
diffusion. The spreading of a stream of blood cells in 
channel flow is characterized by a sub-diffusive behavior 
with exponent 1/3, a phenomenon expected to be generic 
to systems where advected particles undergo short range 
pairwise hydrodynamic interactions or collisions. For 
RBCs, this scaling still holds at significant local concen- 
trations where multibody interactions have to be con- 
sidered in other systems such as rigid beads. We provide 
previously unpublished values of the down-gradient shear 
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FIG. 4: Effective diffusion coefficient A as a function of 
No/h in a flat channel, (a): 2h = 53 um and 2h = 101 um 
(data restricted to initial maximal concentration between 3 
and 16%, resp. 2 and 12 % and No < 3.7 um). Full line shows 
linear fit. (b): 2h — 101 um. Data are extended to initial 
maximal concentration 30% and No < 19.8 um (curve with 
the largest slope in Fig. [3]). Full line shows quadratic fit . 
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FIG. 5: (color online) Effective diffusion coefficient f(yo) = 
SV2Ah/(27R 2 N ) in the 190x99 um 2 channel as a function 
of lateral position yo. f is expected to converge to when 
2/o—^0. Full line shows fit to theory. 



induced diffusivities of RBCs, with a marked difference 
between the diffusivity fi in the direction of shear and fs 
in the vorticity direction. This strong anisotropy should 
be explained by a detailed analysis of RBC collisions at 
the microscopic scale. Our study pertains to dilute to 
semi-dilute suspensions of RBCs for which the connec- 
tive lift flux due to lift of isolated cells and the diffusive 
lift flux can be readily balanced to get the concentra- 
tion profile of the suspension in channel flow. At higher 
hematocrits though, the screening effect on the lift due 
to other cells as well as the local rheology of the concen- 
trated suspension which modifies the flow profile should 
be considered. 
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Supplemental Material for "Lift and 
down- gradient shear-induced diffusion in Red 
Blood Cell suspensions 

In this supplemental material, we derive the advect ion- 
diffusion equation in a 3D Poiseuille flow. The resulting 
simplified equation in the case of a distribution invariant 
in one direction is similar to the one obtained in a simple 
shear flow. We then discuss the existence of self-similar 
solutions for these equations, whose shape strongly de- 
pend on the initial conditions considered. The solution 
relevant to our problem is then detailed. 



ADVECTION-DIFFUSION IN A POISEUILLE 
FLOW 

We consider the stationnary experiment where an ini- 
tial concentration 0(0, z) is continuously injected at the 
inlet x = of a rectangular channel of section 2d x 2h 
in the (y, z) plane. The origin for (y, z) coordinates is 
at the channel center. We neglect diffusion in the flow 
direction compared to advection. The diffusive flux in a 
simple shear flow J = —DV(j) is different whether con- 
centration gradient is in the velocity gradient direction or 
in the vorticity direction. In the first case, the diffusion 
coefficient D reads D = f2R 2 \ A i\(j) an d in the second case 
it reads D = f^R 2 ^^, where R is the typical particle 
size. In particular, the diffusivity is proportional to the 
rate of collisions We assume that the particles ve- 

locity is independent from the concentration and we will 
generally consider that they flow with the same velocity 
u(y, z) as the fluid, which they do not perturb. Expres- 
sions for Stokes flow in rectangular duct can be found in 
Ref. [42]. 

The deformation tensor G is obtained from the velocity 
u(y,z): 



G = 




(4) 



If (xi,xs) are the local coordinates in the (ui,u 3 ) 
frame, the diffusive flux is by definition 



The stationary advection-diffusion equation is: 



u(y,z) 



dx 



dy y dz z 



(6) 



(7) 



where the J components J y and J z are obtained from 
preceding equation by standard variable manipulation: 



J ii 



U(p 2 U<p 



UyU Z ) 



+/2(-£«: 



dz v dy 
2 , 9<P 



UyU z ) 



(8) 



(9) 



dz' x 1 dy UyUz 



)]■ 



If the distribution is initially z-invariant, we expect 
that it will remain so; J z — and §f = then implies: 



»<->^,|«r-vg> + sv4<»rvg) 

(10) 

If the channel is flat, that is d —> oo, the flow is 
parabolic in the z direction, and u y — so with no fur- 
ther assumption, we get: 



u{z)-=R h \u z \-{4>-) 



(11) 



For a more general section, if we consider a narrow 
cloud centered on y$ , we can make the simplification con- 
sisting in considering only the leading order in the devel- 
opment of v and u y around y$. The first one is constant, 
and the second one is constant if yo ^ or of order 1 if 
2/o = 0. Forgetting about this latter very specific case, 
we find: 



where 



du 



m , * -= y, z. Writing T = 

we find that in an orthonormal basis (ui,U2,U3), G is 
transformed into: 



(5) 



Locally, the flow can be described as a simple shear flow 
with shear rate T > 0, where U2 = (1,0,0) is the flow 
direction and U3 = ^(0, u y ,u z ) gives the local velocity 
gradient direction, ui = ^(0, — u z , u y ) gives the vorticity 
direction,. 




u{yo,z)—=R z — , TT^ir)- (I 2 ) 



dx 



Note that u y and u z are also functions of yo . If we 
integrate in the z direction, we get equation: 



where 



Kvo) = R 



dx =XM dy^ 



2 V \J U l+ U z \J U l+ U zJ 

Ju(y ,z)dz 



(13) 



(14) 
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will contribute directly to the effective diffusion coeffi- 
cient measured when considering the widening of a cloud 
of particles. For the flat channel, X(yo) = ^fl a t = fe^W 
and in the general case where thickness 2h would not be 
too large compared to width 2d, lim^^o A (2/0) — ^flat* 



Eq. 



13| is similar to the one that would be obtained 
in a shear chamber with flow in the x direction. In such 
an experiment, the distribution is usually ^-independent 
and its evolution with time can be considered. The shear 
plane is xy and we have : 



d $ d2 , |.| 9 f ± a <l>s 

m =Rm d y ^d y ) 



(15) 



It is equation 13 with x replaced by time t and A 

f2R 2 \il 



SELF-SIMILAR SOLUTIONS 



A slightly more general case than equation 13 is consid- 
ered for further discussion. In case of the usual Brownian 
diffusion, there is no <j) term in the diffusion coefficient. 
In case of 3-body interaction, we have a (j) 2 term. For per- 
fectly spherical particles, 2-body interactions do not lead 
to any diffusion in creeping flows, and 3-body interac- 
tions is the simplest mechanism leading to irreversibility 
and diffusion. If one of these kinds of diffusion is pre- 
dominant, we have equation: 



dx dy^ dy J 



(16) 



where m = 0, 1, 2 and <j) has been rescaled by some typical 
concentration of the problem O , y by A and x by A/^q 1 . 

In Ref. [34], the shear-induced diffusion of platelike 
particles in a flat channel is considered. Two chan- 
nels meet at the inlet of the main channel, one being 
charged in particles with concentration </>o, the other be- 
ing particle-free. Then this Heaviside-like initial distri- 
bution flattens as it is advected along the channel. As 
the data are restricted to the case where concentrations 
near the walls haven't changed significantly (that is to 
say, the observations are made not too far in the chan- 
nel), the authors suggest to modelize their experiment by 
the spreading of an Heaviside distribution, with <j) = <po 
at one side whatever x and (j) = on the other side. In 
other words, they consider the lateral wall is a source of 
particles. 

In that case, it is natural to look for a solution under 
the form (f)(x,y) = ^(77), where 77 = yx~ a . If such a 
solution exists, it will mean that there is a broadening of 
the distribution with a width that grows with x with an 
exponent a. 



Another typical diffusion experiment is the spreading 
of an initial amount of particles initially injected with sta- 
tionary distribution. In that case, if the width increases 
as x a , particles number conservation implies that the am- 
plitude should decrease as x~ a . We shall therefore look 
for solutions under the form (j)(x,y) = x~ a ^(j]). Finally, 
we include both cases by considering <p(x : y) = x~ aUl $(r]), 
with v = 0, 1. 

Equation [16] becomes: 

x- 1 -^(a^ + ar?^ / )+x- a(2+ ^ +m ^ ) (^ m ^ / ) / = 0, (17) 

which yields an equation for ^ as a function of 77 (and 
not x) iff: 



(18) 



in which case there is a possibility for self-similar <\> with 
a typical width scaling as x 

Solutions with reservoirs of particles at one end (y = 0) 
always spreads with exponent 1/2 as in the experiment 
by Rusconi and Stone. On the contrary, the spreading 
of a fixed quantity of particle depends on the diffusion 
process considered: the exponent is 1/2 for Brownian 
diffusion (and the solution of equation 16 is the Gaussian 
profile), but 1/3 for shear-induced diffusion due to pair 
interactions. One sees here that the necessity to have 
particle collision to get diffusion (m > 1) leads to slower 
diffusion (subdiffusion), as neighbors are required and the 
latter are more and more diluted. 3-body interactions 
will lead to an exponent 1/4. 

Thus, the scalings and shape of solutions strongly de- 
pend on initial conditions, a direct consequence of the 
non-linearity of the advection-diffusion equation. 



ANALYTICAL SOLUTION FOR m = 1. 



When v = 1, polynomial solutions of Eq. 16 can be 
found, under the form ^(77) = or ^(77) = ~^rj 2 + 6, 
where b is free. A solution with parabolic profile ^(77) = 



max(0, — ^77 2 



16 



might have 



b) can then be found. Eq. 
other solutions, but numerical resolution with Mathe- 
matica software starting from any reasonable or even 
exotic concentration profile showed convergence towards 
this self-similar parabolic profile, as illustrated in Fig. [6] 
Back to the initial units, the half- width at half- height 
of the concentration profile can be expressed as 



w(x) = w (l + 



9X(yo)N ,1/3 
x) 



4>/2' 



(19) 



where No = j (j){x,y)dy is a conserved quantity and 
wo = w(0) the initial half- width. At long distance from 
origin, w does not depend any more on wq. The exponent 
1/3 as well as the dependency of the effective diffusion 
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FIG. 6: (Color online) Evolution of concentration profile 
along a channel obtained by numerical resolution of equa- 
tion §^ = J^(0§^) with initial condition <j){x — 0,y) = 



dx 



dy ) 



1 + sin(50y)) exp ( - 100?/ 2 



coefficient on the number of particle Nq are direct evi- 
dences for diffusion due to pair interactions. 



Note also that the parabolic profile has steeper edges 
than the Gaussian profile found for Brownian diffusion, 
which is another indication of the necessity to have 
neighbors to diffuse. In the case of 3-body interactions 
[m = 2), a solution ^f(rj) = max(0, \ \fh — r] 2 ) can be 
found, and the edges are even steeper. The exponent is 
1/4 and the effective diffusion coefficient is proportional 
to 7V 2 . 



Note that the relative importance of pair and 3-body 
interactions are controlled by the local concentration, 
while the widening laws depend oniVo, a global variable 
that can hide different situations: a wide and diluted 
cloud or a narrow and concentrated one. 



